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Abstract 

Deformation  mappings  are  considered  that  correspond  to  the  motions  of  lattice  defects,  elastic  stretch  and  rotation  of  the 
lattice,  and  initial  defect  distributions.  Intermediate  (i.e.,  relaxed)  configuration  spaces  associated  with  these  deformation  maps 
are  identified  and  then  classified  from  the  differential-geometric  point  of  view.  A  fundamental  issue  is  the  proper  selection  of 
coordinate  systems  and  metric  tensors  in  these  configurations  when  such  configurations  are  classified  as  anholonomic.  The 
particular  choice  of  a  global,  external  Cartesian  coordinate  system  and  corresponding  covariant  identity  tensor  as  a  metric 
on  an  intermediate  configuration  space  is  shown  to  be  a  constitutive  assumption  often  made  regardless  of  the  existence  of 
geometrically  necessary  crystal  defects  associated  with  the  anholonomicity  (i.e.,  the  non-Euclidean  nature)  of  the  space.  Since 
the  metric  tensor  on  the  anholonomic  configuration  emerges  necessarily  in  the  definitions  of  scalar  products,  certain  transpose 
maps,  tensorial  symmetry  operations,  and  Jacobian  invariants,  its  selection  should  not  be  trivialized.  Several  alternative  (i.e., 
non-Euclidean)  representations  proposed  in  the  literature  for  the  metric  tensor  on  anholonomic  spaces  are  critically  examined. 
©  2003  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  notion  of  a  relaxed  “intermediate”  or  “natu¬ 
ral”  configuration  in  finite  deformation  anelasticity 
or  elastoplasticity,  wherein  each  local  crystal  volume 
element  exhibits  a  stress-free  state,  was  forwarded 
by  many  researchers  in  the  mid-20th  century  [1—8]. 
Such  a  relaxed  configuration  may  correspond  to  the 
intermediate  configuration  arising  from  the  usual 
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multiplicative  decomposition  of  the  deformation  gra¬ 
dient  [9,10],  with  the  relaxed  state  associated  with 
the  local  unloading  of  each  volume  element  from 
its  stressed  state  in  the  current  (i.e.,  deformed  or 
Eulerian)  configuration.  Alternatively,  the  relaxed 
configuration  may  correspond  to  local  unloading  of 
each  volume  element  from  its  reference  (i.e.,  initial  or 
Lagrangian)  state;  such  a  natural  configuration  will 
differ  from  the  initial  configuration  when  the  body 
contains  a  distribution  of  internal  residual  stress  fields 
associated  with  crystal  defects,  for  example  [1 1,7,12]. 

It  is  important  to  note  that  in  most  cases  the 
stress-free  configuration  is  only  locally  coherent  (i.e., 
simply  connected  or  holonomic),  since  each  volume 
element  that  undergoes  relaxation  will  deform  during 
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the  unloading  process  in  a  fashion  that  is  incompat¬ 
ible  with  the  relaxation  deformations  of  its  neigh¬ 
boring  elements.  In  other  words,  a  crystalline  body 
supporting  a  distribution  of  heterogeneous  internal 
stress  fields  associated  with  defects  must  generally 
be  cut  into  multiple  pieces  in  order  to  simultaneously 
relieve  the  stress  in  every  piece.  The  ensemble  of 
pieces  comprises,  in  differential-geometric  terminol¬ 
ogy,  at i  anholonomic  space  [13,14],  or  a  space  that 
does  not  admit  a  homeomorphism  (i.e.,  an  invertible, 
bi-continuous,  one-to-one  mapping)  to  coordinates 
of  a  three-dimensional  Euclidean  space,  the  latter 
being  the  geometric  space  occupied  by  the  coherent 
crystalline  body  prior  to  local  stress  relaxation.  It  is 
also  understood  stresses  are  relaxed  to  zero  only  in 
an  average  sense  in  the  so-called  “stress-free”  con¬ 
figuration.  In  other  words,  when  traction  is  removed 
from  the  external  surface  of  each  local  volume  el¬ 
ement,  self-equilibrating  forces  may  still  be  present 
internally,  within  each  element.  Total  relaxation  of 
all  internal  forces  (e.g.,  atomic  and  sub-atomic  level 
interactions)  can  only  be  achieved  at  observation 
scales  more  far  more  refined  than  those  describable 
by  continuum  elastoplasticity  theories. 

A  well-known  fact  from  differential  geometry  is 
that  for  a  manifold  (i.e.,  configuration  space)  to  ad¬ 
mit  a  single  global  Cartesian  coordinate  frame  whose 
basis  vectors  are  tangent  to  curves  on  the  manifold, 
the  space  must  be  Euclidean  [13,15].  Necessarily  as¬ 
sociated  with  the  Euclidean  property  of  the  space  is 
vanishing  of  the  Riemann-Christoffel  curvature  ten¬ 
sor  constructed  from  partial  derivatives  of  the  metric 
tensor  of  the  coordinate  system,  with  the  metric  tensor 
corresponding  to  the  covariant  identity  tensor  (i.e.,  a 
covariant  Kronecker’s  delta)  when  a  single  Cartesian 
frame  is  selected.  On  the  other  hand,  when  the  space 
is  non-Euclidean  (e.g.,  anholonomic),  prescription  of 
the  Cartesian  metric  tensor  implying  a  single  set  of 
three  orthonormal  basis  vectors  spanning  the  global 
anholonomic  space  may  be  an  inappropriate  assump¬ 
tion  when  such  basis  vectors  are  identified  as  tangents 
to  the  material  manifold  in  the  intermediate  configu¬ 
ration  [16],  In  recognition  of  this  issue,  alternative  co¬ 
variant  defonnation  measures  (e.g.,  elastic,  plastic,  or 
total  strain  tensors)  have  been  implemented  as  metric 
tensors  for  the  purpose  of  lowering  indices  on  con- 
travariant  quantities  referred  to  anholonomic  spaces 
[17-20], 


In  this  work  we  first  review  the  criteria  for  labeling 
a  space  as  either  holonomic  or  anholonomic  and  show 
how  such  spaces  arise  in  finite  elastoplasticity.  We 
next  demonstrate  reasons  why  a  metric  tensor  on  the 
anholonomic  configuration  must  often  be  specified:  to 
define  the  transpose  operation  for  mixed-variant  ten¬ 
sors,  to  conduct  symmetry  operations  on  nominally 
mixed- variant  tensors,  to  define  Jacobian  invariants  of 
deformation  mappings,  and  to  define  scalar  products 
of  vectors  and  contravariant  tensors  referred  to  the  an¬ 
holonomic  space.  Finally,  several  alternative  choices 
for  the  anholonomic  configuration  coordinate  system 
and  associated  metric  tensor  are  presented  and  then 
critically  evaluated  within  the  context  of  finite  elasto¬ 
plasticity. 


2.  Anholonomic  configurations  in  finite 
elastoplasticity 

Consider  the  following  multiplicative  decomposi¬ 
tions  for  the  total  deformation  gradient  F  [21]: 

F  =  FeFp  = 

rpa  rpC  rpP  ^ — V  /i\ 

t  A  ~  r  jx.r  A  —  Jt  aJt  o  A  ,  (l) 

where  the  associated  configurations  and  local  tangent 
spaces  are  illustrated  in  Fig.  1,  and  where  all  indices 


Fig.  1.  Configurations  and  tangent  maps. 
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a,  A,  tx,  a,  =  1,2,3.  We  label  the  reference  config¬ 
uration  (i.e.,  the  initial  state)  as  BK f,  the  current 
configuration  (i.e.,  deformed  state)  as  Bcm,  the  inter¬ 
mediate  configuration  of  finite  elastoplasticity  (i.e., 
the  relaxed  state  achieved  via  locally  unloading  from 
the  current  configuration)  as  BmU  and  the  natural 
configuration  (i.e.,  the  relaxed  state  achieved  via  lo¬ 
cal  unloading  from  the  reference  configuration)  as 
08.  Local  unloading  from  the  reference  configura¬ 
tion  is  associated  with  relaxation  of  internal  stress 
fields  corresponding  to  dislocations  and  other  defects 
present  within  the  crystal(s)  at  the  initial  time,  when 
the  entire  crystalline  body  is  typically  assumed  to  be 
free  of  external  traction.  Local  unloading  from  the 
current  configuration  corresponds  to  simultaneous 
relaxation  of  internal  stress  fields  due  to  dislocations 
and  other  defects  present  at  the  current  time  as  well 
as  any  traction  applied  to  the  external  surface  of  the 
crystalline  body.  As  mentioned  already,  the  partic¬ 
ular  relaxed  configuration  attained  depends  strongly 
upon  the  size  of  each  local  volume  element,  and  the 
so-called  “internal  stress  fields”  that  are  relaxed  here 
may  be  viewed  as  stresses  arising  from  the  traction 
applied  externally  to  the  surface  of  each  local  volume 
element.  The  elastic  and  plastic  tangent  maps  of  mul¬ 
tiplicative  elastoplasticity  are  written  in  Eq.  (1)  as 
Fe  and  Fp,  respectively.  The  tangent  map  Jfo  links 
the  natural  and  reference  configurations,  while  the 
tangent  map  Xf  links  the  natural  and  current 
configurations.  We  remark  that  while  some  authors 
have  multiplicatively  decomposed  our  if  into  elastic 
and  plastic  parts  [22,12],  we  find  the  decomposition 
(1)  more  convenient  in  the  present  setting,  since  in 
our  approach  the  covariant  leg  of  Fp  is  referred  to 
the  Euclidean  space  Ble f.  In  precise  terms,  when  de¬ 
fined  in  terms  of  stress  relaxation,  the  (inverse  of 
the)  elastic  tangent  map  Fe  is  characterized  only  up 
to  an  arbitrary  rotation  tensor  [10].  Some  aspect  of 
the  microstructure  is  needed  to  specify  the  elastic 
(or  plastic)  rotation,  such  as  local  lattice  orientation 
in  classical  crystal  plasticity  theories,  for  example 
(cf.  [12]). 

We  require  in  the  present  work  for  F  to  be  compat¬ 
ible  (i.e.,  integrable,  holonomic,  or  a  true  deformation 
gradient): 


F  = 


fix 

fiX’ 


A  8XA’ 


(2) 


where  x  and  X  are  current  and  reference  coordi¬ 
nates,  respectively,  the  former  assumed  to  be  smooth 
single-valued  functions  of  the  latter  and  time  t.  Fol¬ 
lowing  directly  from  (2),  certain  integrability  condi¬ 
tions  are  automatically  fulfilled: 

fiV  fiV 

A'B  -  8XA8XB  ~  8XB8XA 

=  F\a  ->  F\b  -  F\a  =  2 F\a  b]  =0,  (3 ) 

where,  as  shown,  the  subscripted  comma  denotes  par¬ 
tial  coordinate  differentiation  and  the  bracketed  indices 
are  anti-symmetrized. 

In  contrast  to  F,  the  deformation  maps  Fe,  Fp,  .if , 
and  Cff  o  generally  do  not  fulfill  conditions  analogous 
to  (2)  and  (3);  i.e.,  these  maps  are  generally  not  in¬ 
tegrable.  In  fact,  the  lack  of  their  integrability  is  of¬ 
ten  associated  with  the  presence  of  crystal  defects  (cf. 
[4,5]).  For  example,  consider  the  line  integral  of  the 
differential  vector  element  dx  =  Fp  dX  about  a  closed 
contour  c  in  Blnl: 

/?  =  jj  dx‘  =  j  F(  dXA 

c  C 

=  -  [  f(am  dXA  A  dXB,  (4) 

4  A 

where  we  have  assumed  Cartesian  coordinates 
XA  A  fi,c r  and  invoked  Stokes’  theorem  (cf.  [13,23]) 
to  convert  from  a  line  integral  to  a  surface  integral 
over  the  area  A  enclosed  by  the  loop  C  (the  image 
of  c  in  AVer ),  and  where  the  skew-symmetric  rank  2 
differential  area  element  is  denoted  by  dXA  A  dXB.  In 
Eq.  (4),  b  is  typically  called  the  net  Burgers  vector 
of  all  geometrically  necessary  dislocations  piercing 
the  area  A  at  the  current  time.  Equivalently,  we  may 
write  for  b 


C  C 


=  -  [Fla,l]dxaFdxb  (5) 

J  a 

in  terms  of  the  inverse  of  the  tangent  map  Fe  and  the 
area  element  dxa  A  dxb  in  the  current  configuration. 
We  see  from  (5)  that  the  anholonomicity  of  Bmt  is 
directly  related  to  defects  present  within  the  crystal  at 
the  current  time,  since  configuration  Bml  is  defined  in 
terms  of  local  elastic  unloading  from  Bcm.  It  should 
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be  noted  that  while  Eq.  (4)  defines  b  in  terns  of  dif¬ 
ferentiation  and  line  and  area  elements  referred  to  the 
reference  configuration,  b  of  (4)  is  equivalent  to  b  of 
(5)  and  is  not  indicative  of  internal  stress  fields  due 
to  dislocations  present  at  the  initial  time.  Instead,  we 
must  define  a  different  net  Burgers  vector  6  quantify¬ 
ing  the  anholonomicity  of  configuration  88  and  rep¬ 
resentative  of  the  defect  densities  present  within  the 
crystal  at  the  initial  time  (since  88  is  defined  in  terms 
of  local  elastic  unloading  from  the  initial  state  Bvcl ): 


C  C 


r 


-  [  ^d/Ad/, 

J  a 

j>  &xa=  <j>  jt dxA 

c  C 


(6) 


=  ~  f  A  dXB,  (7) 

J  A 

where  c  is  a  closed  circuit  in  the  natural  configuration 
88.  In  the  interest  of  brevity,  subsequent  equations  will 
focus  upon  the  multiplicative  decomposition  F =FeFp 
and  the  associated  anholonomic  space  Bml;  however, 
our  arguments  will  generally  also  apply  to  the  decom¬ 
position  F  =  Jf  Xfi  1  and  the  corresponding  anholo¬ 
nomic  space  88. 

Basis  vectors  locally  tangent  to  the  Euclidean  spaces 
BTef  and  Bcm,  respectively,  are  written 

S  S 

Ga  ~  8X*’  ga=d^-  (8) 

Corresponding  dual  bases  cotangent  to  Bre{  and  Bcm, 
respectively,  are  written  as  GA  and  ga,  and  are  defined 
such  that  the  following  dual  (i.e.,  scalar)  products  are 
satisfied: 


(GA,GB)  =  %,  (ga,gh)=Sab,  (9) 

where  <5^  and  8ab  are  Kronecker’s  delta  symbols,  i.e. 
3^  =  I  lor  A  =  B.  dAt  =  0  for  A  ^  B,  with  analogous 
relations  for  <Y‘h.  From  (8)  we  see  that  the  basis  vectors 
of  the  Euclidean  spaces  Bre f  and  Bcm  are  holonomic 
basis  vectors,  satisfying  the  integrability  conditions 

[13]  . 


(K» 


The  deformation  gradient  is  then  written  as  follows  in 
terms  of  components  and  basis  vectors: 

r)xa 

F=FaAga®GA  =  —Aga®GA.  (11) 

Similarly,  we  write  for  Fe  and  Fp: 

Fp=/^g*®G^,  Fe  =  Fe“xga  0  gx ,  (12) 

where  the  anholonomic  basis  vectors  g7  and  associ¬ 
ated  covectors  g'  must  obey 


(r,g/»>=^ 


(13) 


in  order  to  ensure  satisfaction  of  the  multiplicative 
decomposition  ( 1 ),  and  where  <5^  =  1  for  a  =  [3  and 
Sp  =  0  for  a  /  / i .  In  general,  g\yji]  0  (compare 
with  Eq.  (10)  for  holonomic  bases),  since  we  cannot 
differentiate  with  respect  to  anholonomic  coordinates 
xx  G  Bint  in  the  conventional  manner  and  since  no  re¬ 
lations  analogous  to  (8)  apply  to  the  ga.  If  we  define 
anholonomic  coordinate  differentiation  as  [13] 


di)=d(lFV-8  =  dAlFe 

8xx  ~  8XA  x  8xa 

then 


(14) 


:  g[M^P«  1  1  =  g[ «,aFep]  £  0 


(15) 


are  additional  conditions  associated  with  the  anholo¬ 
nomicity  of  configuration  Bmt.  On  the  other  hand,  only 
when  Bmt  is  simply  connected  ( and  Euclidean),  mean¬ 
ing  that  the  xx  are  holonomic  coordinates  of  xa  or  XA, 
do  we  have 


F 


e— 1“ 
.a 


8xx  ~  _  8 

fa*’  8“ - 


g[ot,/i]  =  0. 


(16) 


3.  The  metric  tensor  on  an  anholonomic  space 

The  metric  tensor  on  the  anholonomic  space  /)mt  is 
written  here  as  g  =  gapgx  0  gj’‘  and  is  symmetric  by 
definition,  i.e.  gxp  =  gpa.  It  satisfies 

gyp  =  gy  ■  gfl,  (17) 

where  denotes  the  scalar  product  or  inner  product 
of  contra  variant  vectors.  The  scalar  product  of  two 
arbitrary  vectors  a  and  b  on  the  tangent  space  of  B ;nt 
then  becomes,  from  (17), 

a  •  b  =  dx ga  •  Pgp  =  axbp(ga  g p)  =  3*gxpP.  (18) 


G[A,B\  =  0,  g  [a.b]  =  0- 
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Components  of  the  metric  tensor  g  are  used  to  lower 
indices  in  the  conventional  manner,  i.e., 

ax  =  ga/sap,  (19) 


and  are  needed  to  compute  the  transpose  maps  of 
mixed- variant  tensors  (cf.  [15]),  e.g., 

FT  =  =  g^Ff  "gba, 

Ff  =  g^BGAB  =  galtFf  *Gba,  (20) 

where  g  XP  are  components  of  the  inverse  metric  g~ 1 , 
Gab  are  components  of  the  inverse  of  the  metric  ten¬ 
sor  G  on  the  reference  configuration  BTe f,  and  gab  are 
components  of  the  metric  tensor  g  on  Bcm.  Addition¬ 
ally  in  (20)  the  dual  map  operation  corresponding  to 
a  horizontal  exchange  of  indices  is  denoted  by  the  su¬ 
perposed  [24,25].  From  relations  (20),  we  are  in 
position  to  define  symmetric  elastic  and  plastic  strain 
tensors  referred  to  the  current  and  reference  configu¬ 
rations,  respectively,  e.g.  the  covariant  Ce  and  Cp: 


Ceab  =  Fela  gbcFeM  =  C\ 


ha ’ 


C%=FfGcAFvl  =  ClA. 


(21) 


The  metric  tensor  g  is  also  needed  to  symmetrize 
nominally  mixed-variant  tensors  referred  to  configura¬ 
tion  Bmt.  For  example,  consider  the  so-called  “plastic 
velocity  gradient”  Lp: 

Lp  =  FpFp” 1 ,  Zp"  =  FPAFP~  1  ’ ,  (22 ) 

where  the  superposed  dot  denotes  the  material  time 
derivative.  The  transpose  of  the  covariant  version  of 
Lp  is  written 

L%=L%  =  gPxKK~lA  (23) 

from  which  follow  the  covariant  symmetric  rate  of 
plastic  deformation  Dp  and  the  skew-symmetric  “plas¬ 
tic  spin”  Wp: 

2 np  —  fp  _i_  f p  2WP  —  I p  —  Jp 


L%=D%+Kp-  (24) 

Multiplying  (24)  by  through  by  g XP  results  in 
Lpl  =  g  +  g  WPxp  =  Dp ]  +  Wp J,  (25  ) 

where  Dp^  =  gaxDp ^  and  JVpfl  =  gaxWPp  both  depend 
upon  the  metric  g  (from  (23))  and  its  inverse  (from 


(25)).  We  thus  see  that  the  plastic  velocity  gradient 
Lp  can  only  be  decomposed  into  stretch  rate  and  spin 
terms  after  a  metric  tensor  g  has  been  introduced,  as 
noted  by  Maugin  [18]. 

Jacobian  invariants  of  deformation  mappings  Fe  and 
Fp  defining  the  relationships  between  volume  elements 
dr  C  Bcm,  dr  C  Bml,  and  d  V  C  BTe f  also  depend  upon 
g.  For  example,  letting  z’  and  ZA  denote  coordinates 
referred  to  local  Cartesian  frames  on  Bml  and  BK f, 
respectively,  we  write  [8] 


8zx  \ 
~8ZA) 


=  det 


8jF_  dXB\ 
dxP  3  8ZA  ) 


=  det(FpB)^J  det(gxlj)/det(GAlj),  (26) 


where  we  have  used  the  identities  det(gaJ5)  = 
(det(dzx/dxP))2  and  det(GAB)  =  (det(8ZA/8XB))2  (cf. 
[15]).  Analogously  for  Fe  we  write 


dr 


dr 


8za 


8zx 


Je  =  —  =  det  —  =  det  — ■  Fc„  — 


8za 


dxP 


dxh  P  dz« 


=  det  (F^ )  ^det(gab)/det(gap),  (27 ) 

with  za  denoting  Cartesian  coordinates  on  BCUI  such 
that  det(grfli)  =  (det(&fl/5x*))2.  We  can  alternatively 
write  (26)  and  (27)  as  [26] 

JP  =  l^BCFP>pBFp( 

=  \{JWG)e^eABCFp)FiFPL  (28) 

r  =  l^SahcF'lF^F* 

=  6  ( Jg/g  )erJ,l/'  eabcFc“  f),  F&z  ,  (29) 

where  we  have  used  the  shorthand  notation  g  = 
det(gu),  g  =  det (gap),  and  G  =  dct( G/Ui ).  Permu¬ 
tation  tensors  in  Eqs.  (28)  and  (29)  are  defined  by 

^  \  G C  '  ,  C,y[jy  =  \/ gCyjly  ,  C>ah(  =  \/geab(y 

and  =  ~\fgexPx,  with  eABC,  expx,  eabc ,  and  e *PX 

standard  permutation  symbols  each  referred  to  the 
configuration  indicated  by  its  indices.  Consider  for 
example  the  incompatibility  described  by  Eq.  (5), 
which  can  be  rewritten  as 

=  -  [ F\~^  dxa  Adxb  =  ^  [  sabcFl~i;nc  da 
J  a  J  a 


sabcF 


e— 1“ 
.a,b 


nc  da 


Ae  nc  da, 


(30) 
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where  the  second  rank  area  element  has  been  con¬ 
verted  to  an  axial  covector  n  da  via  dx"  A  dxh  = 
sabcnc  da,  and  where  the  two-point  dislocation  density 
tensor  Ae  =  —sahcFeN  .  Mapping  the  integrand  in 
(30)  to  the  intermediate  configuration  via  Nanson’s 
formula  (cf.  [27] )  gives 

bx  =  f  Ae  ncda  =  I 'jeFe~1'  Ae  rip  da 
J  a  J  a 

=  I  Ae'  rip  da,  (31) 

J  a 

where  Ae‘  =  JeFe~1'  Ae  is  the  intermediate  configu¬ 
ration  dislocation  density  tensor,  which  depends  upon 
Je  and  hence,  from  Eq.  (29),  g.  Should  we  opt  to  allow 
the  (scalar)  residual  free  energy  i//<lN  of  the  crystal  as¬ 
sociated  with  geometrically  necessary  dislocations  to 
depend  quadratically  upon  a  dislocation  density  mea¬ 
sure  referred  to  Bmt  such  as  Ae  (cf.  [27-29]),  we  then 
write,  for  example, 

i£GN  =  i/fr2 Ae :  Ae  =  \nl2Ae’f  ,  (32) 

where  fi  and  /  are  an  elastic  shear  modulus  and  char¬ 
acteristic  length,  respectively.  Thus,  the  scalar  prod¬ 
uct  in  Eq.  (32)  is  another  instance  in  finite  crystalline 
elastoplasticity  where  g  is  required. 


4.  Metric  tensors  and  anholonomic  basis  vectors: 
possible  choices 

We  now  discuss  some  conceivable  choices  for  the 
metric  tensor  g  and  associated  basis  vectors  ga.  The 
simplest  and  by  far  most  prevalent  option  from 
the  literature  is  to  impose  (cf.  [7,8,22,30]) 

(jctfl  =  ga  '  g/j  =  brjfj,  (33) 

with  djj;  the  covariant  Kronecker’s  delta,  i.e.  8ap  =  1 
for  ot  =  fi  and  =  0  I'or  a  /  fi.  Eq.  (33 )  indicates  that 
contravariant  vectors  defined  on  Bmt  are  referred  to  a 
single  global  Cartesian  frame  (i.e.,  orthonormal  ga), 
or,  equivalently,  a  parallel  (i.e.,  identical)  Cartesian 
frame  attached  to  each  local  relaxed  volume  element 
du  C  2?int-  Such  a  choice  implies  either  that  (i)  Bmt  is  a 
Euclidean  space  or  that  (ii)  the  ga  are  not  actually  tan¬ 
gent  to  any  coordinate  lines  scribed  on  the  base  man¬ 
ifold  ^jnt  (since  such  coordinates  do  not  exist  in  the 


discontinuous  anholonomic  space)  but  instead  corre¬ 
spond  to  some  external  reference  frame.  As  discussed 
already  in  Section  2  of  the  present  work,  statement  (i) 
is  ruled  out  in  the  presence  of  crystal  defects  such  as 
dislocations  that  render  Bmi  anholonomic.  This  leaves 
statement  (ii),  which  is  difficult  to  interpret  geomet¬ 
rically:  from  its  standpoint,  vectors  in  configuration 
2?jnt  are  not  referred  to  the  actual  tangent  spaces  of  a 
material  manifold  in  Bmt  but  are  instead  referred  to 
the  fixed  external  Cartesian  frame(s).  The  intermedi¬ 
ate  metric  8xp  (or  equivalently,  the  external  Cartesian 
frame)  is  introduced  rather  artificially  in  such  theo¬ 
ries  as  an  additional  constitutive  entity,  accompanying 
variables  Fe  and  Fp.  We  also  mention  the  related  work 
of  Simo  [31],  wherein  the  intermediate  configuration 
metric  is  given  by  g.,p  =  <5^ Gab 8Bp,  implying  a  Eu¬ 
clidean  metric  structure  for  Bml  when  Bref  is  Euclidean, 
and  relying  upon  the  external  two-point  construct  8Ay. 
While  we  are  unable  to  rigorously  rule  out  the  choice 

(33)  on  any  fundamental  grounds,  we  emphasize  next 
several  alternatives  that  have  been  proposed  in  the 
literature  that,  in  contrast  to  Eq.  (33),  do  not  make 
the  embedding  of  a  non-Euclidean  space  Bmt  within 
a  global  Cartesian  space  equipped  with  a  Euclidean 
metric  tensor  <5ajg.  (It  should  be  noted,  however,  that 
our  derivations  (26)  and  (27),  wherein  Cartesian  za 
are  assigned  to  each  volume  element  du  C  /imt,  do 
rely  on  such  an  assumption,  at  least  locally. ) 

One  such  alternative  is  to  specify  the  metric  ten¬ 
sor  g  as  the  covariant  elastic  deformation  tensor  Ce 
[17,18,20],  i.e. 

9^  =  Clp  =  Fcl_gabFcl  =  Fe°x ga  ■  F^gb  =  ga  ■  g/j, 

(34) 

implying  that  the  anholonomic  basis  vectors  are  de¬ 
fined  by  ga  =  Fcrj  g„ .  Analogously,  one  could  imple¬ 
ment  the  covariant  plastic  defonnation  tensor  Cp  as  a 
metric  [32] 

9.P  =  C%=F^AGabfY“ 

=  F*~ lA Ga  ■  Fp~iBGb  =  ga  -gp,  (35 ) 

n_lA 

meaning  that  g7  =  F,,  GA.  (It  should  be  noted  that 
Miehe  [32]  suggested  several  alternative  metric  ten¬ 
sors  on  each  configuration,  in  addition  to  (35 ). )  In  Eqs. 

(34)  and  (35)  we  have  used,  respectively,  the  standard 
Euclidean  relations  gab=ga'gb  and  Gab=Ga-Gb-  The 
get  in  (34)  or  (35)  are  still  not  tangent  to  any  global 
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coordinate  curves  x7  (since  such  coordinates  are  pre¬ 
cluded  by  the  anholonomicity  of  /}]nl ),  but  they  do  ex¬ 
ist  in  a  one-to-one  manner  with  the  set  gfl(x)  and  the 
map  Fe  or  with  the  set  G.ffX)  and  the  map  Fp  1 ,  re¬ 
spectively.  It  follows  that  Ce  and  Cp  are  well-defined 
geometric  quantities  referred  to  the  space  J3jnt,  in  con¬ 
trast  to  the  somewhat  artificial,  external  object  dy/i  of 
(33)  that  cannot  be  derived  from  the  elastic  or  plastic 
tangent  maps  and  basis  vectors  on  BTe f  or  Bcm.  Con¬ 
sider  now  the  rank  4  Riemann-Christoffel  curvature 
tensor  formed  from  Ce  on  the  anholonomic  space  Bmt, 
with  components  [13,14] 


pa  — 

Kpxs- 


r)fx  rtf 7 

01  &t>  _  01  xP 

dx%  dxs 


r.p 


IB  W 


rc  4-  ?i?c  r* 

1  Se 1  x 


(36) 


where  the  coefficients  of  the  metric  connection  are 
found  by  2 C^f  ^  =  Cc/M  +  Cesp  x  -  C%a  and  compo¬ 
nents  of  the  anholonomic  object  by 
The  conditions  R°p  s  =  0  and  k?s  =  0  hold  identically 
only  when  Fcrj  =xaa,  i.e.,  only  when  Bmi  is  homeomor- 
phic  to  the  Euclidean  space  Bcm  and  x7  are  holonomic 
coordinates  for  which  we  can  define  partial  differen¬ 
tiation  as  usual.  Under  such  conditions,  KJ;./A  is  the 
pull-back  of  the  curvature  tensor  formed  from  the  met¬ 
ric  tensor  gab,  which  itself  vanishes  since  Bcm  is  Eu¬ 
clidean.  Analogous  statements  regarding  the  curvature 
can  be  made  for  Cp,  i.e.  the  curvature  tensor  derived 
from  Cp  vanishes  identically  only  when  F\  =x7/l  and 
B i„t  is  Euclidean. 

A  third  alternative  methodology  [33]  involves  in¬ 
troduction  of  the  new  multiplicative  decomposition 

F  =  FeFp,  FaA  =  Fc'‘HFpl,A,  (37) 


where  both  legs  of  the  new  plastic  deformation  tensor 
Fp  are  referred  to  /iref  and  where  Fe  is  a  two-point 
tensor  between  tangent  spaces  of  BIe f  and  Bcm,  i.e. 

Fe  =  Fe^ga  <g>  Gb ,  Fp  =  Fp"aGh  ®  Ga.  (38) 


Following  Le  and  Stumpf  [33],  we  introduce  the 
two-point  transformation  matrix  H,  satisfying 

B 


such  that  Eqs.  ( 1 )  and  (37)  are  satisfied  simultaneously 


TP  a  /-’e"  r-r  p  pe"  rj  —  1  t_tCC  rpP  pe''  tt’P 

^  A  ~  r  A  ~  r  .B^i  a  H  £t  A  —  t  Bt  A. 


(40) 


The  anholonomic  basis  vectors  and  covectors  are 
found  in  terms  of  FI  as,  respectively. 


ga=H~l  Ga 


r  =  o7aga. 


(41) 


so  that  the  dual  product  obeys  the  standard  relations 

(£«,&)  =  {h:/ga,h^gb)  =  h~v‘  h^(ga,gb) 

=  H-iAH<ibSb  =  SP,  (42) 

leading  then  to  the  representation  of  the  metric  tensor 
on  Bmt: 

(jocll  =  ga  •  g IS  =  H  y  GA  ■  H  p  G B 

=  H-'AGAliH-f.  (43) 

Eq.  (41)  implies  that  the  elastic  deformation  maps  Fe 

and  Fe  are  the  same  (two-point)  tensor,  although  each 
is  referred  to  different  coordinate  bases: 

Fe  =Fe; ga  ®  g7  =  (KHf  )ga  ®  (H«b Gb) 

=  Fe“.Bga®GB  =  Fe.  (44) 


Analogously,  for  the  plastic  tangent  mappings  we  have 

Fp  =  F? ga  ®  G^  =  (Fp>rc)(77-|,iG«)  ®  G^ 

=  F'p(4Gb®Ga  =Fp.  (45) 

From  the  anholonomicity  conditions  (15)  we  see  that 

§[«./(]  =H\*,p]GA  +hT)  ga,P]  =H~{i I^Ga 


+  (46) 

Upon  choosing  a  Cartesian  reference  coordinate  sys¬ 
tem  for  the  Euclidean  space  Bref,  we  obtain  GA  B  =  0, 
and  Eq.  (46)  becomes 

g[*/n  =  nGA  *  0  -  vfn  *  «,  (47) 

meaning  that  H  (or  its  inverse)  is  generally  not  deriv¬ 
able  as  the  gradient  of  a  vector-valued  function  when 
Bmt  is  anholonomic.  On  the  other  hand,  only  when  the 
intermediate  configuration  is  holonomic  and  g[a,fl  =0 
can  we  write 


dy* 


dXA 


-1 
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where  y  1 = y  7(X'\  t )  are  new  coordinate  functions  on 
Bmt,  not  necessarily  identical  to  xa.  Upon  invocation 
of  the  holonomicity  relations  (16)  we  obtain 

,e»  dy  a  dxa  dy  a 


FeA=Fe, 


dXA  dx*  dXA  ’ 


_  dXB  «  _  SXB  dx* 

A  dy*  A  dy a  dXA  ’ 


(49) 


from  which,  for  the  particular  choice  of  x*  =  y  “,  re¬ 
sults  in  FeA  =  FaA  and  FVA  =  SBA.  Returning  now  to 
the  anholonomic  case  and  substituting  (39)  into  (4), 
we  see  that  the  incompatibility  (i.e.,  closure  failure  or 
net  Burgers  vector  in  configuration  /}ml)  depends  si- 

„  B 

multaneously  upon  gradients  of  F  A  and  H*c\ 

b«  =  -  f  (HFp)a[A  B]  dXA  A  dXB 
Ja 


=  -  [ (F^  B]H«c+FpC[AHyB])dXA AdXB,  (50) 
J A 

implying  that  knowledge  of  components  and  gradi- 

zv  B 

ents  of  Fa  alone  is  not  sufficient  to  characterize  the 
incompatibility  of  Eq.  (4)  (i.e.,  one  needs  either  FVA 

(Eq.  (4)),  Fc-V‘  (Eq.  (5)),  or  both  of  FPA  and H*B  (Eq. 
(50)).  Thus,  while  attractive  at  first  glance  since  no  an¬ 
holonomic  basis  vectors  directly  enter  the  multiplica¬ 
tive  decomposition  (37),  the  methodology  outlined  in 
Eqs.  (37)-(50)  can  only  be  used  within  the  context 
of  classical  higher-order  gradient,  dislocation-based 
crystalline  elastoplasticity,  wherein  the  incompatibil¬ 
ity  is  referred  to  the  anholonomic  intermediate  con¬ 
figuration  Rjnt,  when  the  additional  tensor  H  and  its 
inverse  are  available  for  mapping  to  and  from  that  con¬ 
figuration.  However,  the  method  should  definitely  not 
be  disregarded,  as  Le  and  Stumpf  [33]  successfully 
applied  this  description  to  define  a  different  Burgers 
vector  B  (i.e.,  an  incompatibility)  referred  to  the  holo- 
nomic  reference  configuration,  in  terms  of  the  lack  of 
integrability  of  Fe,  i.e. 

BA  =  f  Fe~ lA  dxa  =  -  j F*r*  d/Adx*.  (51) 

C 

It  should  be  noted  that  in  this  context,  components  of 
FeA  are  not  associated  with  loading  from  a  locally  un¬ 
stressed  and  anholonomic  state  to  a  stressed  state,  but 
rather  describe,  at  each  point  in  the  crystal,  the  stretch 
and  rotation  of  a  triad  of  lattice  director  vectors  occur¬ 
ring  between  the  initial  configuration  and  the  current 


configuration,  with  Fp  =Fe  1 F  any  remaining  “inelas¬ 
tic”  deformation  referred  to  the  reference  state. 


5.  Discussion 

Additional  critical  comments  regarding  the  par¬ 
ticular  choices  of  intermediate  configuration  metric 
tensor  listed  in  Eqs.  (33),  (34),  (35),  and  (43)  are 
now  in  order.  Consider  first  the  designation  given 
by  (33):  gy  =  by,.  This  prescription  is  by  far  the 
most  common  in  finite  deformation  plasticity  theories 
from  the  literature,  as  it  is  made  implicitly  whenever 
Cartesian  coordinates  for  the  intermediate  configura¬ 
tion  Bmt  are  invoked,  be  it  for  the  purpose  of  defining 
differential  volume  change  [7,8],  calculating  a  norm 
or  scalar  product  of  contravariant  vectors  or  tensors 
[27-30],  or  decomposing  the  plastic  velocity  gradi¬ 
ent  Lp  into  distinct  deformation  rate  and  spin  terms 
(cf.  [34]).  We  remark  that  Regueiro  et  al.  [29]  ac¬ 
knowledged  the  non-Euclidean  character  of  Bmt  but 
gave  no  explicit  alternative  to  the  Cartesian  metric 
tensor  for  use  on  Blnl.  As  stated  already,  choice  (33) 
cannot  be  ruled  out  unequivocally  via  mathemati¬ 
cal  or  physical  arguments.  However,  since  the  unit 
metric  tensor  (i.e.,  covariant  Kronecker’s  delta  dy) 
is  not  considered  a  well-defined  geometric  variable 
on  a  non-Euclidean  space  (cf.  [16]),  and  since  Bmt 
is  clearly  non-Euclidean  (i.e.,  not  homeomorphic  to 
three-dimensional  Euclidean  space)  under  conditions 
of  heterogeneous,  incompatible  plastic  deformation, 
the  metric  dy  of  (33)  must  be  viewed  as  an  addi¬ 
tional,  stationary  constitutive  variable  introduced  by 
the  modeler,  accompanying  any  constitutive  relations 
specifying  the  time  evolution  of  components  of  Fe 
and  Fp.  In  contrast,  the  choices  given  in  (34),  (35), 
and  (43) — wherein  the  intermediate  metric  tensors 
are  defined  directly  in  terms  of  elastic  or  inelastic 
tangent  mappings — may  be  viewed  as  favorable  al¬ 
ternatives  to  (33)  since  the  former  do  not  require 
the  modeler  to  prescribe  an  intermediate  metric  via 
introduction  of  an  additional  constitutive  variable. 
Maugin  [19]  voiced  similar  opinions  regarding  the 
choice  of  metric  tensor  on  our  natural  configuration 
(Fig.  1 ),  cautioning  against  usage  of  a  Cartesian  met¬ 
ric  on  a  non-Euclidean  space.  Instead,  Maugin  [19] 
elected  to  employ  the  quantity  XarigahXfh A  as  a  met¬ 
ric  for  lowering  contravariant  indices  on  the  generally 
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anholonomic  space  'A.  It  should  also  be  noted  that 
while  the  covariant  object  cannot  be  constructed 
directly  from  anholonomic  mappings  Fe  or  Fp,  the 
mixed-variant  unit  tensor  (i.e.,  true  Kronecker’s 
delta)  can  be  introduced  explicitly  via  the  inverse 

operation,  i.e.  Safj  =  Fc-r Fc“fj  or  dajs  =  FpAFpfA . 

Eq.  (34),  prescribing  gap  =  Ceap  =  F^gahFcp  as  a 
possible  alternative  to  (33),  has  been  proposed  ex¬ 
plicitly  by  at  least  several  others  [17,18,20],  and  is 
analogous  to  the  aforementioned  approach  favored  by 
Maugin  [19]  when  Jf  is  replaced  by  Fe.  The  quantity 
is  considered  to  be  a  well-defined  geometric  entity 
with  indices  referred  to  the  intermediate  configuration 
B[nt,  since  its  constituents  Fca  and  gab  are  by  defini¬ 
tion  single-valued  functions  of  current  coordinates  xa . 
(This  is  in  contrast  to  the  external  construct  Sxr  which 
must  be  introduced  by  an  additional  constitutive  as¬ 
sumption.  )  Moreover,  the  Riemann-Christoffel  curva¬ 
ture  tensor  RpxS  and  the  anholonomic  object  Kj,y  of 
Eq.  (36)  both  vanish  identically  only  when  Fe_1  is  an 
integrable  function  of  x  and  Smt  is  Euclidean,  mean¬ 
ing  that  assignment  of  Cyp  as  components  of  a  metric 
tensor  agrees  fully  with  the  non-Euclidean  nature  of 
the  generally  anholonomic  space  Sjnt.  Consider  also 
the  implications  of  choice  (34)  in  a  thermodynamic 
assumption  for  stored  elastic  energy  such  as  Eq.  (32): 
the  free  energy  associated  with  geometrically  neces¬ 
sary  dislocations,  written  herein  as  t//GN,  would  depend 
explicitly  upon  the  elastic  strain  metric  C®^.  Such  an 
assumption  is  in  fact  supported  by  previous  physical 
[35]  and  numerical  [36]  experiments  on  various  duc¬ 
tile  pure  metals  and  their  alloys,  wherein  amplification 
of  internal  energy  of  dislocation  arrays  at  flexing  grain 
and  sub-grain  boundaries  in  the  presence  of  applied 
loads  was  discovered.  In  light  of  the  above  reasons, 
we  endorse  here  usage  of  C®^(  as  a  viable  alternative  to 
£>o:/j  as  a  metric  tensor  on  the  globally  non-Euclidean, 
anholonomic  space  /imt. 

Prescription  of  the  plastic  strain  tensor  = 
Fp~l  GahF^p  '  as  a  metric  gxp  was  suggested  in  Eq. 
(35),  following  previous  work  by  Miehe  [32],  The 
quantity  CPp  is  considered  to  be  a  well-defined  ge¬ 
ometric  entity  with  indices  referred  to  the  interme- 
diate  configuration,  since  its  constituents  Fpa  and 
Gab  are  by  definition  single-valued  functions  of  refer¬ 
ence  coordinates  XA  and  time  t.  This  approach  appears 


especially  attractive  since  for  the  case  of  no  plastic 
deformation  (and  hence  no  defects  or  anholonomic- 
ity),  we  have  Fp~x  =  SAy  and  gap  =  SAaGAB <5^,  thus 
reducing  to  the  prescription  of  Simo  [31].  However, 
from  the  perspective  of  Eq.  (32),  this  method  presents 
some  difficulties,  since  the  plastic  strain  CPp  would 
enter  the  free  energy  expression  i//GN,  and  since  we 
do  not  view  the  plastic  strain  in  most  cases  as  a  useful 
state  quantity  or  accurate  measure  of  stored  elastic 
energy.  The  latter  point  is  of  no  concern  in  local 
theories  (e.g.  Miehe  [32])  wherein  higher-order  de¬ 
formation  gradients  are  not  explicitly  considered  and 
energetic  quantities  such  as  i//GK  of  (32)  are  not  com¬ 
puted.  But  in  nonlocal  thermomechanical  theories 
framed  in  the  relaxed  intermediate  configuration  (cf. 
[27-29])  wherein  relations  such  as  (32)  are  invoked, 
then  we  regard  the  choice  of  intermediate  metric  C®^ 
(34)  as  more  physically  plausible  than  the  choice 
=  C%  (35). 

Finally,  consider  the  aforementioned  approach 
of  Le  and  Stumpf  [33],  with  the  multiplicative  de¬ 
composition  in  component  form  given  by  Eq.  (37): 

*  a  *  B 

FaA  =  FCHF'A .  As  discussed  already,  this  decom¬ 
position  enjoys  the  simplicity  of  not  requiring  any 
anholonomic  basis  vectors  or  corresponding  interme¬ 
diate  metric  tensor.  Moreover,  this  method  offers  an 
alternative  viewpoint  of  kinematics  of  elastic  defor¬ 
mation  in  terms  of  the  total  deformation  (between  Sre f 
and  BCUI )  of  a  triad  of  lattice  director  vectors  assigned 
to  each  material  point.  However,  as  we  have  shown 
in  Eq.  (50),  quantification  of  the  anholonomicity  of 
Sint  (i.e.,  the  lack  of  integrability  of  FPA)  requires 
knowledge  of  the  anholonomic  mapping  H,  which 
in  turn  implies  the  intermediate  configuration  metric 

iA  i  B 

g,jjt  =H~l  GABH  p  as  noted  in  Eq.  (43).  (It  should 
be  mentioned  that  the  metric  H~  GABH  ^  is  gener¬ 
ally  non-Euclidean,  so  this  choice  does  not  contradict 
the  non-Euclidean  nature  of  the  anholonomic  space 
.Sint.)  From  the  perspective  of  Eq.  (50),  since  HJH 

r,  B 

and  its  gradient  are  required,  in  addition  to  FA  and 
its  gradient,  this  method  is  more  complicated  than 
the  other  choices  (34)  and  (35),  which  require  only 
higher  gradients  of  either  F\  (Eq.  (4))  or  Fe~v  (Eq. 
(5)).  This  aspect  leads  us  to  favor  g =  C®^  (34)  over 
gxp=H-lAGABH-lB  (43). 
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6.  Concluding  remarks 

In  higher-order  gradient  plasticity  theories  couched 
explicitly  upon  the  mechanics  of  incompatible  defects, 
special  care  should  be  taken  when  prescribing  metric 
tensors  for  the  anholonomic  (e.g.,  natural  or  interme¬ 
diate)  configurations  that  are  inherent  to  the  kinematic 
description,  since  various  scalar  products,  transpose 
maps,  tensorial  symmetry  operations,  and  Jacobian 
invariants  depend  explicitly  upon  the  chosen 
anholonomic  basis  vectors  and  their  inner  product.  We 
discuss  several  alternative  choices  to  the  popular  and 
apparently  contradictory  (but  not  incorrect)  choice 
of  a  Euclidean  metric  tensor  (and  accompanying 
Cartesian  basis  vectors)  attached  to  an  anholonomic 
(and  hence  non-Euclidean)  configuration  space.  En¬ 
dorsed  here  as  an  intermediate  configuration  metric 
is  the  pull-back  of  the  current  configuration  metric 
tensor  by  the  elastic  deformation  mapping  (i.e.,  the 
covariant  elastic  strain  referred  to  the  intermediate 
configuration). 
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